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Abstract 

We report the discovery of a new mechanism of spontaneous generation of a magnetic flux in 
a superconductor cooled through T^. The sign of the spontaneous flux changes randomly from 
one cooldown to the next, and follows a Gaussian distribution. The width of the distribution 
increases with the size of the temperature gradient in the sample. Our observations appear incon- 
sistent with the well known mechanisms of flux generation. The dependence on the temperature 
gradient suggests that the flux may be generated through an instability of the thermoelectric 
superconducting-normal quasiparticle counterflow. 
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With the exception of ferromagnets, a spontaneous appearance of a magnetic field in a 



physical system is a highly unexpectec 
superconductors cooled through T, 



lenomenon. Yet, such phenomenon was observed in 
. In one case, the spontaneous magnetic 




field appeared as a consequence of the d-wave symmetry of the order parameter of HTSC 



_cooling the superconductor 
. Here, we report a new. 



. In another case, the spontaneous field was generated 
through Tc under non equilibrium thermal conditions 0, P. 
unexpected appearance of a spontaneous field, which occurs in a superconductor cooled in 
the presence of a thermal gradient. 

Although the effect described below appears completely unrelated, the original motivation 
for this experiment followed our previous work on the Kibble-Zurek cosmological scenario j3, 
0]. One of the key assumptions of this scenario is that the temperature within the sample 
is uniform. The limit on the size of VT, the temperature gradient across the sample, set by 
Kibble and Voloviky], is that VT < Tci/C,. Here is the transition temperature, i = 
and ^ are the reduced temperature and coherence length respectively at the temperature T at 
which fluctuations of the order parameter return to thermal equilibrium ^. Our experiment 
was designed to see what happens to the formation of topological defects once this criterion 
is not satisfied. 

The experimental setup is the same as described in Ref. 4, with the exception of a 
non-uniform heating, generating intentional temperature gradients in the sample. Briefly, 
our samples were 300 nm thick c-axis oriented YBa2Cu307 films with Tc ~ 90 K, grown 
on a SrTiOs substrate. The samples were placed atop the sensing coil of a HTSC SQUID 
magnetometer. In our arrangement the SQUID remains at a temperature of 77 K, and is not 
affected by the temperature of the sample which can be heated and cooled independently. 
The ai,n i= heated above u=ing a hght sou.ce and cools by exchanging heat with .ts 
environment. The light source is a pulsed YAG laser [10[. Single pulses (FWHM ~ 10 ns) 
were used to heat the film. The laser pulse passes through the substrate and illuminates 
non-uniformly a selected area of the film. At a laser wavelength of 1.06 /xm, the SrTiOs 
substrate is transparent and practically all the light is absorbed in the film. Hence, only the 
film heats up, while the substrate remains near the base temperature of 77 K. The 1 mm 
thick substrate has a heat capacity about 10'^ larger than that of the film. The heat from 
the film escapes into the substrate, which acts as a heat sink. This small thermal mass of 
the film allows us to achieve cooling rates in excess of 10*^ K/sec. The cooling rate at 



2 



can be varied by changing the amount of energy dehvered by the laser pulse. The system is 
carefully shielded from the earth's magnetic field, with a residual field of less than 50 fiG. 
Additional small coil adjacent to the sample was used to test the field dependence of the 
results, at fields ranging from less than 50 fiG up to 60 mG. 

Non-uniform illumination was generated either by using a non-uniform light beam, or by 
covering some part of the sample. An example of one such arrangement is shown in the 
inset of figure |21 Here, the strongly illuminated area is a stripe across the film. In another 
configuration, the perimeter of the film was masked, while an area of 4mm in diameter in 
the center was exposed to the beam. Qualitatively, the results presented here do not depend 
on the exact illumination profile. Under such non-uniform illumination, the film cools down 
in a two stage process. In the first stage, the heat deposited by the laser pulse in the film 
is dumped into the SrTiOs substrate on a time scale of a fis. As a result, the temperature 
of the part of the substrate closest to the illuminated area increases by up to 5K above 
77K, depending on the laser energy. In the second stage, heat is transferred from the hotter 
part of the substrate to its cold parts on a time scale of several tens of ms. During all this 
time temperature gradients are present across the sample. The relatively slow time scale on 
which the substrate cools is due to its thermal mass, which is much larger than that of the 
film. 

Previously, under homogeneous illumination (VT ~ IK/cm), we have observed the gen- 
eration of spontaneous fiux during a rapid quench of a superconducting filmj^. The flux 
appeared faster than the temporal resolution of our SQUID, which is ~ 10/is. In the follow- 
ing, we refer to this signal as the "fast" signal. The polarity of the flux from one quench to 
the next was random, following approximately a Gaussian distribution centered at zero. The 
width of this distribution increased weakly with the quench rate, a result which is broadly 
consistent with the Kibble-Zurek scenario. 

Under non-homogeneous illumination, we estimate that VT increased to about 300 K/cm 
for the largest pulse energy used. This is still less than the limit set by the homogeneous 
criterionyi of 10^ K/cm. Under these conditions, the "fast" signal showed no appreciable 
change. However, in addition to the "fast" signal, an unexpected, much larger signal has 
appeared after a relatively long delay of 1-10 ms (see flgure^). This signal was completely 
absent during measurements using a homogeneous illumination. We point out that the time 
it takes to cool below Tc is on the order of 1 fis. Consequently, the "slow" signal appears 
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while the film is already in the superconducting state. 

The polarity of the non-homogeneous, "slow" signal was also random from one quench 
to the next. Similarly to the "fast" signal, the amount of flux generated in a given quench 
followed a Gaussian distribution centered at zero. This is shown in Figure 2. However, the 
amount of flux associated with the "slow" signal is larger than that of the "fast" signal by 
an order of magnitude (see also figure E)) . 

After analyzing data aquired using different pulse energies, we found that the amount of 
spontaneous flux, characterized by the distribution width, increases with the pulse energy. 
This contrasts the results found under the conditions of uniform illumination (there the 
distribution width decreased with increasing pulse energy). This is clearly seen in figure 01 
in which the signal dependence on pulse energy is shown. Note that increasing the pulse 
energy also increases the thermal gradients generated across the film. 

Finally, measurements were repeated under different external magnetic fields ranging 
from less than 50 /iG up to 60 mG. As figure El clearly shows, the results do not depend on 
the external field. 

The results at non-homogeneous conditions point towards two important conclusions. 
First, as already noted above, increasing the temperature gradients across the film by 2 
orders of magnitude (from 1 K/cm up to 300 K/cm) does not change the "fast" signal. 
Therefore we conclude that the homogeneous approximation 9] indeed holds, at least for 
thermal gradients up to ~ 10^ K/cm. Second, the dependence of the "slow" signal on pulse 
energy and the long time scale clearly imply that it originates from another mechanism, 
rather than the Kibble-Zurek scenario. In the following, we discuss several other mechanisms 
which may generate magnetic flux, and examine their possible relevance to our observations. 

predicts a conversion of thermal energy into magnetic 
field fluctuations while the sample is in the critical region near T^. In our experiment, the 
sample passes through this region in less than 1 /is, while the slow signal develops on a time 
scale 3 to 4 orders of magnitude slower, 1-10 ms. So, this scenario does not fit with our 
observations. 

Another possibility is a change in the spatial distribution of residual magnetic flux inside 
the film. Re-arrangement of magnetic flux lines can happen during partial illumination of 
the samples. Magnetic flux can move in or out of the heated part of the film, changing 
the magnetic flux distribution. Re-distribution of magnetic flux can then change the actual 
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amount of flux coupled to the SQUID, even though the net change is zero. We investigated 

this mechanism in separate measurements done at the university of Konstanz, Germany using 

I I 

a magneto-optic system capable of sub ns resolution [12]. We found that re-distribution of 
flux takes place within several ns, which again is inconsistent with the time scale of our slow 
signal. In addition, re-distribution of flux should depend on the ambient fleld, which is not 
borne by our data. 



Several theory papers 



14| proposed that flux can be generated by an instability of a 



propagating normal-superconducting phase boundary front, which indeed is present in our 
samples as the fllm cools after a non-homogeneous heating pulse. If this mechanism is viable, 
it should act during less than 1 /is after the heating pulse, since at later times the entire fllm 
cools back into the superconducting state and the front disappears. Again, this is 3 order 
of magnitude faster than the time after which the slow signal is observec . 

Spontaneous flux can be formed at large angle grain boundaries jl|, 15] as a consequence 
of the d-wave symmetry of the order parameter. Our samples are epitaxial thin fllms, in 
which large angle grain boundaries are absent. Hence, this mechanism can not explain the 
origin of our signal. 

One clue as to the origin of the effect comes from the observation that the temperature 
gradients across the sample relax on the same time scale as the time over which the slow 
signal develops. Therefore it is natural to associate it with some thermoelectric effect. This 
association would also be consistent with the size of the effect increasing with the energy 
deposited in the fllm. Thermo-electric effects (the Seebeck effect or the Nernst effect jlfil|) 
can generate flux lines as a result of superconducting currents in the fllm. 

In superconductors, the Nernst effect is a result of the motion of flux lines along the 
thermal gradient. Clearly, this effect depends on the ambient magnetic fleld. Since we see 
no such dependence, we conclude that the Nernst effect does not explain our measurements. 

Regarding the Seebeck effect in a superconductor, thermal gradients produce a counter- 
flow of normal quasiparticles and Cooper pairs. The net electric current is zero However, 
as noted by Ginzburg^^, in some cases such thermo-electric currents can generate magnetic 
flux. One example is , in which the supercurrent 

and the normal current are not co-linear and form a current loop. This happens if the 
Seebeck coefficient is anisotropic and the direction of the thermal gradient is not parallel to 
one of the superconductor's symmetry axes. Then, the superconducting countercurrent does 
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not exactly cancel the normal current at every point of the film, hence generating a non- 
zero magnetic flux. Measurements done by Subramaniam et al. [l9| show that for untwinned 
YBCO crystals, thermoelectric properties are indeed anisotropic. However, our films are 
twinned, so there is no anisotropy between the a and b directions (parallel to the surface 
of the film). Under a temperature gradient of 300 K/cm, we estimate the thermo-electric 
current / ~ 5 x lO^^A. This estimate is based on the measured thermal coefficients (l9|. 

If the spontaneous flux was generated via a linear thermoelectric effect, we would expect 
the polarity of the flux generated to be the same in each measurement, since the temperature 
gradient in the sample is nominally the same. Since the polarity of the measured flux is 
random in each measurement, this suggests that perhaps an instability occurs. 

One well known example is the plasma "two stream instability" '2d|. In a plasma, the 
Lorentz force between the two opposing electron beams vanishes only if the currents cancel 
exactly everywhere. With spatial current fluctuations, the cancellation does not hold, re- 
sulting in a net repulsive force between the currents. This in turn leads to further separation 
of the currents, creating a current loop and a magnetic flux. This situation is illustrated 
in figure lU The sense of the current in the loop (clockwise or counterclockwise) is ran- 
dom, having been determined by the initial fluctuation which separates the current and 
countercurrent. This is in line with our data. In superconductors, the counterfiow consists 
of a normal current opposed by a supercurrent. A high frequency plasma instability in a 
superconductor was predicteded by Kempa et al. |2l|. Another variant of the "two stream 
instability" was proposed by Bliokh and Shapiro [22], who showed that in the framework of 
the two-fluid model 2^, a uniform superconducting-normal quasiparticle countercurrent is 
unstable with respect to spatial fluctuations. The analysis reveals that fluctuations in the 
current density can induce a low frequency instability and generate a magnetic field. De- 
noting the velocities and densities of the normal and superconducting components by Vn^s , 
Un^s, and taking n^K = n„y„ and + n„ = n, the growth rate of the unstable mode has 
the formQ: 

u; = ^^ (1) 

UnUs V 

where v is the electron relaxation time and k is the wavevector of the unstable mode. 
To obtain a numerical estimate, we determine Vg from = en^K, using the thermoelectric 
current estimated above, / ~ 5 x lO^^A and the sample cross section. The wavevector k is 



determined by taking the size of a fluctuation as A, the penetration depth. This choice fixes 
k — 27r/A. For other quantities in (1), we used the two fluid model expressions, a charge 
density n of 10^^ holes/cm^, and u ~ 10^^ Hz. Assuming the current flows at the final stage 
of the experiment near the boundary of the sample, we find u ~ lO+'^Hz- 10^'^Hz, which is 
consistent with the measured experimental growth rate (ranging between 10 Hz - lO^Hz). 
Hence, this scenario gives a possible explanation for the origin of the measured spontaneous 
flux. 

In conclusion, we have discovered a new mechanism of spontaneous flux generation in 
a superconductor quenched through Tc in the presence of a temperature gradient. One 

mechanism which may be responsible for this new effect is an instability of the thermoelectric 
current distribution. 
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FIG. 1: Typical SQUID signals showing the fast (top trace) and slow(bottom trace) formation of 
spontaneous flux (note the different scale of the horizontal axes.) Arrows show the time at which 
the laser pulse was applied. 

FIG. 2: Typical distribution of spontaneous flux under non-homogeneous illumination. Solid black 
bars show the noise distribution, while the dashed curve shows a Gaussian fit to the signal distri- 
bution. The inset shows typical non-homogeneous illumination profile. 

FIG. 3: Signal distribution as a function of pulse energy, showing the difference between homoge- 
neous and non-homogeneous illumination. Also shown are measurements done at several different 
external magnetic fields. The error bars are statistical. 
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FIG. 4: A schematic picture of the current loop formed by the super and normal thermoelectric 
currents, separated as a result of the instability. 
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